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ABSTRACT 


Tv/o problems  are  considered^  both  of  which  arc  characterized  by  essentially 
the  same  geometry.  In  the  first  problem,  the  properties  of  surface  waves  propagating 
in  a  slab  of  uniaxially  anisotropic  dielectric  are  studied*  The  slab  is  of  finite  thick¬ 
ness  and  the  relative  dielectric  constants  of  the  dielectric  tensor  are  all  positive.  Two 
cases  are  studied  simultaneously,  one  being  the  case  where  the  optical  axis  is  chosen 
parallel  to  the  slab,  and  the  other  being  the  case  where  the  optical  axis  is  chosen  per¬ 
pendicular  to  the  slab. 

Only  £-mode  surface  waves  are  studied  since  these  have  the  properties  of  ex¬ 
traordinary  type  waves.  The  H-mode  surface  waves  (ordinary  surface  waves)  have  the 
same  properties  as  surface  waves  propagating  in  an  isotropic  slab.  Hence,  the  results 
for  the  latter  problem  are  well  known. 

A  transverse  resonance  procedure  is  utilized  in  the  determination  of  the  surface 
wave  resonances.  It  is  found  that  at  any  frequency  the  resonances  are  finite  in  number. 
In  addition  it  is  found  that  surface  waves  will  propagate  only  if  certain  conditions  are 
satisfied  by  the  principal  dielectric  constants  of  the  diagonal  tensor. 

In  the  second  probelm  the  electromagnetic  fields  of  a  uniaxially  anisotropic 
plasma  slab  in  the  presence  of  a  magnetic  line  source  are  studied.  The  geometry  is 
the  same  as  for  the  first  problem,  except  that  the  optical  axis  is  chosen  parallel  to  the 
slab  and  perpendicular  to  the  direction  of  the  line  source.  The  line  source  is  located 
either  outside  the  slab  or  at  a  conducting  plane  (1.  e.  .  an  infinite  slot)  covered  by  the 
slab. 

Under  the  conditions  described,  it  is  found  that  the  structure  can  support  an 
infinite  number  of  discrete  £-mode  surface  waves  for  any  frequency  less  than  the 
plasma  frequency.  The  resultant  surface  wave  field  can  be  expressed  as  a  convergent 
series  obtained  by  evaluating  the  residue  terms  which  arise  in  the  steepest  descent 
representation.  Furthermore  it  is  found  that  all  the  surface  waves  are  forward  waves 
in  so  far  as  propagation  adong  the  slab  is  concerned* 

In  addition  to  the  surface  waves,  an  infinite  number  of  nonspectral  leaky  wave 
poles  is  found  to  exist  for  any  frequency  greater  th^h  plasma  frequency.  However, 
only  a  finite  number  of  the  leaky  wave  poles  contribute  residue  terms  for  any  angle  of 
observation  less  than  90^,  where  the  angle  of  observation  is  measured  from  the  normal 
to  the  slab. 
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CHAPTER  I 
INTRODUCTION 

The  problems  to  be  considered  in  subsequent  chapters  will  involve  uniaxially 

anisotropic  dielectric  media.  Media  possessing  such  properties  can  be  constructed 

by  implanting  dielectric  obstacles  in  a  suitable  dielectric  binder.  Such  obstacles 

would  have  dimensions  which  are  small  in  comparison  to  the  wavelength  of  the 

radiation  and  would  be  arranged  in  a  geometric  pattern  which  would  give  rise  to  the 

uniaxial  properties.  If  the  obstacles  possess  finite  conductivities^  the  resulting 

artificial  dielectric  structure  would  be  anisotropic  both  electrically  and  magnetically. 

This  possibility  will  not  be  considered  here.  The  properties  of  artificial  anisotropic 

dielectrics  as  well  as  the  techniques  for  constructing  them  have  been  considered  by 
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various  authors,  and  the  results  are  available  in  the  literature  *  ’  *  . 

The  characteristics  of  surface  waves  guided  by  isotropic  dielectric  slabs  have 
been  studied  In  great  detail  and  the  results  are  well  known An  extension  to  this 
problem  consists  of  the  study  of  such  waves  guided  by  a  uniaxially  anisotropic  di¬ 
electric  slab.  This  latter  problem  will  be  considered  in  Chapter  III  and  It  will  be 
found  that  the  surface  waves  which  propagate  at  a  given  frequency  are  finite  in  number. 
In  addition,  it  will  be  shown  that  the  relative  dielectric  constants  of  the  dielectric 
tensor  must  satisfy  certain  Inequalities  if  surface  waves  are  to  propagate. 

Another  medium  which  possesses  the  uniaxial  characteristics  described  above 
is  a  plasma  In  the  presence  of  an  infinite#  d.  c.  magnetic  field.  If  a  slab  of  such  a 
plasma,  with  the  optical  axis  oriented  parallel  to  the  slab,  la  excited  by  a  magnetic 
line  source  some  very  interesting  phenomena  are  observed.  These  will  be  discussed 
in  detail  in  Chapter  IV.  However,  it  is  interesting  to  mention  some  of  the  results 
here,  and  by  way  of  contrast,  to  mention  some  of  the  results  obtained  by  Tamir^^* 
in  the  analysis  of  the  source- excited#  isotropic  plasma  slab. 

In  the  case  of  the  uniaxially  anisotropic  slab  it  will  be  shown,  for  the  E-mode 
case#  that  an  infinite  number  of  surface  waves  can  be  supported  by  the  slab  at  any 
frequency  below  the  plasma  frequency.  In  addition,  the  waves  are  all  forward  waves. 
In  contrast,  Tamir*s  results  show  that  an  isotropic  plasma  slab  can  support  up  to 
four  forward  or  backward  surface  waves,  and  these  carry  power  in  opposite  directions 
in  the  plasma  and  air  regions.  Other  contrasting  results  are  also  found  to  exist. 

Above  the  plasma  frequency,  in  the  anisotropic  case.  It  is  found  that  a  descrete 
infinity  of  nonspectral  leaky  wave  poles  exist#  only  a  finite  number  of  which  contribute 
residue  terms  to  the  field.  In  Tamir's  problem  the  nonspectral  leaky  waves  do  not 
occur.  Of  course,  the  continuous  spectrum  associated  with  open  structures  occurs 
in  both  the  isotropic  as  well  as  the  anisotropic  case. 
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CHAPTER  II 

FORMULATION  OF  THE  GENERAL  EQUATIONS 

The  repreientation  of  electromagnetic  fields  in  cylindrical  regions  is  greatly 
simplified  by  expressing  the  fields  in  terms  of  their  longitudinal  and  transverse 
components  with  respect  to  the  longitudinal  coordinate  z.  The  transverse  fields  can 
then  be  represented  in  terms  of  a  complete  set  of  vector  modes  characteristic  of  the 
transverse  field  distribution!  and  the  longitudinal  fields  can  be  derived  directly  from 
the  transverse  fields.  In  this  way  the  problem  is  reduced  to  one  of  solving  a  set  of 
one-dimensional  transmission  line  equations.  The  techniques  are  equally  valid  for 
isotropic!  as  well  as  uniaxlally  anisotropic  regions. 

The  formulation  described  above  will  be  outlined  for  a  region  filled  with  a 
medium  of  permeability  and  a  relative  dyadic  permitlvity  given  by 

e  s 2-0  io  ®1 2^0  *1-0-0 

for  the  case  where  the  optical  axis  coincides  with  the  x  axis  (o  a  x)!  and 

®  o— o  ^  “^l^oJ^o  ^  ®2— o— o^  (2-2) 

for  the  case  where  the  optical  axis  coincides  with  the  z  axis  (o  a  z)«  The  inhomo¬ 
geneous  Maxwell  equations  then  have  the  form 

V  X  E  =  -  H-M  (2-3) 

and 

y  X  H  =  JuuCq  €  •  5.  £  (2-4) 

where  and  ^  «e  the  electric  field!  the  magnetic  field!  the  magnetic  source 

current  density  and  the  electric  source  current  density!  respectively.  A  time 
dependence  e^^^  is  assumed. 

The  Maxwell  equations  can  be  put  into  the  equivalent  transverse  (to  z)  form 

9  19 

by  Using  the  standard  techniques  outlined  in  the  literature.  *  It  is  pointed  out  that 
the  results  will  be  given  for  the  two  dimensional  case  (i.  e.  ^/dy  »  0)1  since  the 
problems  considered  subsequently  fall  into  this  category.  For  the  transverse  fields 
and  E  ^  one  readily  obtains  the  following  equations: 
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’’*]■ 


(H.  X  a  )  + 


X  z 

—  t  -'O 


(2~5) 


where 


=  ju)^ 


f 


t 


(2 


+  Zq 


.2  2 

*"o  'Vo' 


(2-6) 


and 


(o  a  x) 
(o  a  z) 


£t 


'lioio  ^*2X0^0  ' 


The  equivalent  transverse  source  currents  and  are  given  by 

V.  X  (z^  J  J 


4  =M *' 

-t  -t  jwe^Cj 


and 


A 

J  =  J 
-t  -t 


o  *1 

V.  X  (a  M  ) 
t  — o  z' 

i<»Un 


(2-7) 


(2-8) 


The  longitudinal  field  components  can  then  be  expressed  in  terms  of  the  transverse 
field  components*  Thus^ 


and 


H  =-r^ -  rv*.(z  xEJ-M  1, 

a  juJMo  L  *  “ 

E  =T— i - U.  .  (H.  X  a  )  -  J  1. 

a  JiDe^ei  L  *  “*  “J 


(2-9) 


(2-10) 


For  the  problems  to  be  considered  the  transverse  region  is  unbounded  in  the  trans** 
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verse  direction,  so  that  a  radiation  condition  must  be  satisfied  at  infinity.  In  the  z 
direction  the  boundary  conditions  will  be  specified  for  the  particular  geometries 
considered. 

Since  the  region  under  consideration  is  uniform  in  the  z  direction,  the  trans¬ 
verse  components  of  the  electric  and  magnetic  fields  can  be  expressed  in  the  form 
V(z,  e^  (x,  and  I(z,  h  (x,  ^),  where  ^  is  the  continuous  (real)  modal  wave- 
number  in  the  transverse  (x)  direction.  Also,  the  homogeneous  nature  of  the  trans¬ 
verse  cross  section  permits  the  decomposition  of  the  vector  eigenfunction  into  E  and 
H  modes  relative  to  z.  This  is  the  same  as  the  representation  for  the  isotropic  case. 
Thus,  the  transverse  fields  and  source  terms  are  expressed  as  follows: 

OD  09 

E^(x,z)=J  V'(z,?)  e'(x,?)  d?  +  J  V'(z,?)  A'(x,(;)d?.  (2-lla) 


Hj(x,  z)  =  J  I'(z.?)h'(x,J:)  d?+  J  I'(z,c)h'(x,?)d?. 


(2-llb) 


?j(x,  z)  =J  i'(z,  e;)e'(x,  ?)d?  +  J  l'(z,  e)e'(x,  5)d?  , 


(2-llc) 


Mj(x,  z)  =  J  v'(z,  ?)  h/(x,  ?)d5  +  J  v*(z,«:)h'  (x,;)d;«  (2-lld) 

ii  “  £  O  *  —  ^ 

where  the  single  and  double  primes  denote  E  and  H  modes. 

Upon  substituting  equations  (2-11)  into  (2-5)  and  (2-6),  and  remembering  the 
facts  that  for  E  modes  e^  =  M  modes  h  =  h 

e*  are  linearly  independent  vector  functions,  one  obtains  the  transmission  line 
equations  and  associated  eigenvalue  problems.  Thus,  for  the  (o  a  x)  case,  the 
separation  procedure  yields  for  a  particular  E  mode 


dz 

i)=j 

(z,?)  +  v'(z,i;). 

(2-12a) 

i)=j 

k'  VMz,5)  +  1'  (x.c). 

Z 

(2-12b) 

- 1 

^  z'  --t  =  . 

(2-12c) 

^  z  V 

V  o«i  ?  ,  z  , 
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and 


(2-13a) 

VjVj.h'U.?)  =  0. 

(2-13b) 

The  E  modea  are  ueually  referred  to  aa  extraordinary  type  modea. 
same  procedure  for  a  particular  H  mode  one  obtains 

By  applying  the 

=  J  k\Z'I'  (x,5)  + y'  (a,?). 

(2-14a) 

■  *^*-i*-*^-^  =  J  V'  (a,  5)  +  1'  (a.  ?). 

(2-14b) 

-  /  2  2  r,»  O 

a=V''o«l-?  '  ^  "5^ 

(2-14c) 

and 

(2-l5a) 

O 

II 

> 

(2-15b) 

The  H  modes  are  referred  to  as  ordinary  type  modes* 

The  same  procedure  outlined  aboret  applied  in  the  (o  a  s)  case  yields  similar 
results.  Thus  for  £  modes  (extraordinary  type  modes)  one  obtains 

(a,  ?)  +v'  (a.?) 

(2-l6a) 

5)  +  i'  (*.?) 

(2-l6b) 

a  V  €£  V  '2  ^  *  *®o®l 

(2-l6c) 

and 

7j7j  .  e'  (x,5)  =  -  ?^a'  (x,  ?), 

(2 -17  a) 

(2-17b) 

Fop  a  particular  H  mode  (ordinary  type  mode)  the  results  are 


(2-18a) 
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=Jii\  Y'  V'  (*.§)  +!'(*,?) 

(2.18b) 

!  in - T  '"‘^o 

5  •  ^  ■"77“ 

(2-18c) 

(2-19a) 

.  e'  (x,5)  s  0. 

(2-19b) 

It  has  been  mentioned  that  the  problems  to  be  considered  subsequently  are 
charactized  by  unbounded  cross  sections  in  the  transverse  direction*  Hence*  equation 
(2-13a)9  (2-15a)*  (2-17a)  and  (2-19a)  must  all  satisfy  the  radiation  condition  at  infinity* 
Appropriate  solutions  for  these  equations  would  then  be  of  the  form 

*  (x,  g)  =  — —  (2-20) 

for  £  modes*  and 

h(x,e)  =— i-  (2-21) 

for  H  modes* 

The  vector  functions  £  and  h  are  chosen  to  satisfy  the  usual  orthogonality 
conditions^*  so  that  the  various  current*  voltage  and  source  amplitude  functions  are 
given  by  the  following  equations: 


V(z*  5)  =  |^E^(x*  z)  .  £*(x*  ?)  dx  * 

(2-22a) 

I(*.?)  =1^  Hj(x,  *)  .  b*(x,  ?)  dx, 

(2-22b) 

v(e, Aj(x,  *)  .  h*(x, ?)  dx. 

(2-22c) 

l(*i?)  =  J  ^(*.  *)  •  «*<*•?) 

(2-226) 

where’^  denotes  complex  conjugate* 

Although  the  preceding  formulation  has  been  outlined  for  two  dimensional 
uniaxlally  anisotropic  regions*  the  problems  to  be  considered  will  also  Involve  isotropic 
regions*  Hence*  the  modal  formulation  for  such  regions^  will  also  be  utilised  in  sub¬ 
sequent  work* 
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CHAPTER  III 

SURFACE  WAVES  GUIDED  BY  A  UNIAXIALLY 
ANISOTROPIC  DIELECTRIC  SLAB 

In  this  chapter  the  properties  of  surface  waves  guided  by  a  uniaxlally  anisotric 
dielectric  slab  will  be  studied,  with  the  analysis  being  limited  to  the  source  free  case. 
The  propagation  of  surface  waves  in  a  slab  of  dielectric  (isotropic  or  anisotropic) 
depends  upon  the  phenomenon  of  total  reflection  at  the  inner  boundaries  of  the  slab.  It 
is  well  known  that  when  total  reflection  takes  place  at  a  dielectric  interface  the  propaga¬ 
tion  constant  In  the  direction  perpendicular  to  the  interface  and  outside  the  dielectric 
becomes  purely  imaginary.  The  fields  outside  the  dielectric  then  decay  exponentially 
with  distance  from  the  interface. 

The  geometry  of  the  problem  to  be  considered  is  illustrated  in  Fig.  3.1.  The 
slab  consists  of  an  anisotropic  dielectric  with  dielectric  tensor  given  by  equation  (2-1) 
or  (2-2),  and  with  and  ^nd  positive.  The  medium  which  surrounds  the  slab 

on  either  side  is  isotropic  with  a  relative  dielectric  constant  ip.  Because  of  the 
obvious  symmetry  of  the  problem 


O  ■" — o 


Fig.  3. 1:  Anisotropic  dielectric  slab  and  equivalent  transmission  line, 
two  bisections  will  be  considered,  a  short  circuit  or  electric  wall  (s  b),  and  an  open 
circuit  or  magnetic  wall  (o  b).  The  analysis  will  be  limited  to  the  study  of  E  modes 
since  the  properties  of  H  modes  (ordinary  modes)  in  the  slab  are  the  same  as  those 
of  waves  propagating  in  an  isotropic  region. 


1.  The  Transverse  Resonance  Solution 
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It  will  be  convenient  to  anticipate  the  fact  that  the  surface  wave  resonances  are 

finite  in  number  at  any  given  frequency.  Hence,  the  subscript  m  will  be  introduced  on 

all  subsequent  results  to  denote  a  particular  surface  wave.  The  wavenumbers  inside 

and  outside  the  slab  in  the  x  direction  are  denoted,  respectively, by  Hr  and  k 

xm  xm 

These  must  be  equal  for  the  boundary  conditions  to  be  matched  at  the  interface.  The 
characteristic  impedance  of  the  equivalent  transmission  line  is  given  by 


Z  = 


m 


zm 


(o  a  x) 


(3-la) 


for  I  z  I  <  dt  and 


(o  a  z) 


k 

Z  _ _ L 

om  oue 


(3-lb) 


(3-2) 


for  I  z|>  d.  And  are  the  wave  numbers  in  the  z  direction,  inside  and  outside 

the  slab  respectively. 

The  values  of  k  »  h  i  And  k  for  the  surface  waves  which  will  propagate 
xm  zm  zm  * 

are  obtained  by  solving  simultaneously  the  transverse  resonance  equation  and  the 
dispersion  relation.  This  latter  relation  is  obtained  from  equations  (2-12c)  and  (2-l6c) 
for  the  region  |z|<  d.  For  the  present  problem,  the  dispersion  relation  within  the 
slab  is  written  as  follows: 


*^0  ®1  ®2 


r  2  2 

’‘xni®2  ®1 

’‘xm®!  ^’‘zL®2 


(o  a  x). 


(o  a  z). 


(3-3a) 

(3-3b) 


For  the  region  |  z  |  >  d. 


the  dispersion  relation  is  given  by 


1  2  2  ^  u  2 

h  e  =  K  +  h 
o  ^  xm  zm 


(3-4) 


The  transverse  resonance  equation  for  either  bisection  is  given  by 


Z  (d)  +  Z  (d)  =  0, 

4-  “• 

Where  Z  (d)  is  the  impedance  looking  to  the  left  at  z  =  d  and  Z  (d)  is  the  impedance 
looking  to  the  right  at  z  =  d.  Thus,  for  short  circuit  or  opA  circuit  bisections  taken 
at  z  ::  0  in  Fig.  3. 1,  one  obtains  the  following  equations: 
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0  b  case 


o  b  case 


tanK__d+  -5J2.  =0,  (oax)  (3-5a)  cot  k  d+ =  o,  (oax)  (3-6a) 

62  ®  6  2  ® 

tanK^^d  =n,  (oa*)  (3.5b)  cotK,„d+ =  0,  (oaz)  (3-6b) 


n 


zm 


A  necessary  condition  for  the  propagation  of  surface  waves  within  the  slab  is 
that  the  propagation  constant  in  the  z  direction  in  the  outer  medium  be  purely 
imaginary,  that  is 


k  =  -j  k 
zm  ]  zm 


With  this  last  result,  equations  (3-3)  and  (3-4)  can  be  combined  to  yield 
62(61  -  e)  ei+Km|^®2  ' 

for  the  (oax)  case»  and 

k^  61(62  -  e)  =»‘,Le2 

for  the  (o  a  z)  case.  Now,  Introducing  the  change  of  variables 


(3-7) 


(3-8) 


P  *  H_d 
zm 


(3-9) 


and 


‘'zmh  ' 


equations  (3-7)  and  (3-8)  both  take  the  form 


(3-10) 


2^2  2 
p  +  q  =  a 


(3-11) 


where 
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(k  d)^ (e.  -  e)  .  (oax), 

(k  d)  (c2  “  ®)  »  *  *)• 

o  *^2  ^ 


(3.12a) 

(3-12b) 


Under  the  same  change  of  variables  (3-9)  and  (3-10)»  equations  (3-5)  and  (3-6)  become 


p  tan  p  ,  (s  b), 


(3-l3a) 


p  cot  p  »  (o  b). 


(3.13b) 


Note  that  equation  (3-13)  is  valid  for  both  the  (o  a  x)  case  and  the  (o  a  s)  case« 

Equation  (3-11)  is  the  equation  of  a  circle  and  it  is  to  be  solved  simultaneously 
with  equation  (3-13).  The  graphical  technique  for  the  solution  Is  illustrated  in  Fig.  3.2. 

The  cutoff  frequencies  for  the  various  surface  waves  occur  when 
approaches  k^e»  ^  which  case  it  follows  from  equations  (3-4)  and 


Fig«  3.  2:  Illustration  of  the  graphical  technique  for  the  transverse  resonance  solution. 
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{3«10)  that  q  =  0,  To  satisfy  these  conditions  the  relations 

tan  p  =  0  cot  p  =  0 

must  hold  in  the  (s  b)  and  (ob)  case  respectively.  The  last  equations  are  satisfied  when 


rr  d=a  =  22JL 


(3-14) 


where  m  =  0,  1,  2,  •  •  •  Note  that  m  is  even  for  the  (sb)  case  and  odd  for  the  (ob)  case. 
The  cutoff  frequencies  and  wavelengths  can  be  obtained  fromequation8(3- 12)  and  (3-14). 
When  the  optical  axis  coincides  with  the  x  axis  one  obtains 


cm 


4d 


and 


m  c- 


c) 


(3-15) 


(3-16) 


where  c  is  the  speed  of  light  in  free  space.  When  the  optical  axis  coincides  with  the  z 
axis  the  results  are 

_  /  •’n 

(3-17) 


cm 


me 

4d 


and 


cm 


^  4d  jTl 
mV  c" 


(Eg  -  0 


(3-18) 


It  should  be  pointed  out  that  the  lowest  cutoff  frequency  is  zero  In  the  case  of  the  short 
circuit  bisection. 

Some  important  facts  can  now  be  obtained  by  interpreting  some  of  the  above 
results.  Reference  to  equation  (3-12)  shows  that  when  the  optical  axis  coincides  with 
the  X  axis  no  surface  waves  can  propagate  if  <  e*  Likewise,  when  the  dielectric 
is  anisotropic  in  the  z  direction,  no  surface  waves  will  propagate  if  €2  £  c* 


2.  The  Dispersion  Curves 

An  overall  picture  of  the  behavior  of  the  surface  waves  can  be  obtained  from 
an  examination  of  the  dispersion  curves  given  by  vs.  Vd.  X  ^  i«  the  wavelength 
in  the  x  direction.  In  order  to  obtain  the  dispersion  curves  it  will  be  necessary  to 
obtain  an  expression  for  ^/X^^  i®  terms  of  p  and  q.  This  expression,  together  with 
a  knowledge  of  the  behavior  of  the  roots  of  the  transverse  resonance  equation  (see 
Fig.  3.2),  enables  one  to  plot  Vx^  vs.  ^/d. 

Equations  (3-3)  and  (3-4)  can  be  written,  with  the  aid  of  (3-10),  as  follows: 


|o  a  x)  case 
‘2*  " 


(o  a  z)  case 


""  t 

«2 


(3-19b) 
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q2  cj  =  -{V)2  c^c  +  (r  ^d)2  (3.20a)  =  -k^  c^r  +  (X^^d)^  e^.  (3-20b‘ 

Now,  combining  equations  (3-19)  and  (3-20)  solving  for  dividing  by  and  using 

equations  (3-11)  and  (3-12),  one  obtains 

(o  a  x)  case  (o  a  z)  case 


(3-21a) 


t-f- 


2  .  2 
p+  62  q 

'2^  2 
p  +  q 


(3-21b) 


Equations  (3-21)  can  now  be  used  to  obtain  the  qualitative  curves  shown  in  Fig.  3.  3. 

It  is  seen,  from  Fig.  3.2,  that  for  a  infinite,  q  becomes  infinite  and  p  approaches  some 
multiple  of  ^/2,  for  a  particular  surface  wave  resonance.  Hence  all  the  curves 
intersect 


Fig.  3.  3:  Dispersion  curves  for  the  surface  waves. 

the  ^/X  axis  at^TT,  where  g,  is  defined  in  Fig.  3.  3.  The  lowest  surface  wave 
resonance  (m  =  0)  has  an  asymptotic  value  of  *^/Xjj  -  JT"  a  approaches  zero.  This 
fact  can  be  obtained  from  equations  (3-13a)  and  (3-21a).  The  higher  order  resonances 
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have  curves  which  end  on  the  line  ^ ■- fT  at  values  of  Vd  given  by  equations  (3-16) 
and  (3-18). 

The  slope  of  the  curves  is  Fig.  3.  3  is  always  negaUve  and  consequently  the 
surface  waves  are  all  forward  waves  in  the  x  direction.  It  is  expected,  then,  that 
the  net  transfer  of  power  is  in  the  positive  x  direction.  This  is  actually  the  case  since 
the  only  real  component  of  the  Poyntlng  vector  is  P  =  -  E  H  *  (see  equations  (3-26) 
and  (3-27)  in  section  3).  ^ 

3.  The  Spatial  Variation  of  the  Fields 

It  is  a  simple  matter  to  derive  the  field  expressions  for  the  surface  waves. 
These  can  be  expressed  directly  in  terms  of  the  current  function  which  satisfies  the 
equivalent  transmission  line  in  Fig.  3. 1.  It  will  only  be  necessary  to  obtain  an 
expression  for  the  y- component  of  the  magnetic  field  H  ,  the  electric  field  being  given 
by 


1 

(o  a 

x) 

1  1 

(3-22a) 

eoe2 

—  # 

=■-  i 

,  ja  1  <  d, 

1 

(o  a 

*) 

(3-22b) 

1 

aH 

(o  a 

x) 

(3-23a) 

f 

~  4 

1 

(o  a 

>1  ■  1  <  d. 

(3-23b) 

i*«o®2 

» 

1 

t  1 

l.l 

>  d, 

(3-24) 

l‘! 

1 

.1 

>  d. 

(3-25) 

• 

‘1 

If  the  amplitude  of  the  current  function  at  s  =  d  is  specified  as  and  if  one 
keeps  in  mind  the  continuity  requirement  for  and  at  ttie  slab  interface,  then, 
with  tile  aid  of  equations  (3-22)  to  (3-25).  the  following  field  distributions  are  obtained: 
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a  b  case 

(a)  0  <  z  <  d 


T  M  .  Bin  K  2 

^  M  ’‘am  ™ 

*  ^  757  j  ’ 


-U  X  ‘=°*  "am* 

”y^7;=*  T^TirrsT 

TZit 


^d  "xm  ^  -Jx  *  "*m* 


^*"72^r  “'o^i 


e  '^'xm 


C®*  X-rr.^ 
zm 


where 


(b)  z  >  d 


®i  ■  •{ 


ej  ,  (o  a  x) 
.  (o  a  z) 


"am  "^""am*^ 
"  ’/Z7"  ^'"  ®o«j  ‘'®*  "am 


ei  =  < 


€2  t  (o  a  x) 
gj  ,  (o  a  z)  , 


^  -Jk  *  •  (z  -  d) 

-j-  e  •^^xm  e  |  zm|  » 


^d  -J  K  *  "  I  ( z  -  d) 

H - ^  e  ^  ’^xm  e  |  «ml  • 


E  -  A.  *”■■■■ 

'  ^  "'o' 


-ix  *  -  ^  I  (z  -  d) 

e  ‘'^xm  e  zm|  # 


o  b  case 
(a)  0  <  z  <  d 


E  - ^ 

^  JW 

I, 


"zm 


j  ttl  € 


-Jh  X  ‘=‘’*  "zm 
«  ’"‘xm 


o^J 


tin  K,_J 
zm 


a  •  1  »  H 

/I  -  W  *  *ni 

H  =-2-  e  •^^xm  -rj— - 3- 

y  -‘“"zm** 


E  -  .  ^  **xm  - 

a  ■  lUeoCi  “^“^zm 


•*“  x,_z 

zm 

T  • 


(3- 26a) 


(3-26b) 


(3-27a) 
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Equations  (3-26)  and  (3-27)  show  that  the  spatial  variation  of  the  fields  is 
sinusoidal  within  the  slab  and  exponentially  decaying  outside  the  slab.  Furthermore, 

It  can  be  seen  from  Fig.  3.2  that  as  frequency  increases  for  a  given  surface  wave 
the  wave  becomes  more  tightly  bound  to  the  slab.  Fig.  3.  4  shows  the  spatial  varlaHon 
with  B  of  for  the  m=.0  surface  wave.  The  field  distribution  is  shown  for  a  low 
frequency  as  well  as  for  a  high  frequency. 
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CHAPTER  IV 


THE  ELECTROMAGNETIC  FIELD  OF  A  UNIAXIALLY  ANISOTROPIC 
PLASMA  SLAB  EXCITED  BY  A  MAGNETIC  LINE  SOURCE 


The  problem  to  be  coneldered  next  consists  of  an  analysis  of  the  properties 
of  the  electromagnetic  field  of  a  unlaxlally  anisotropic  plasma  slab  excited  by  a 
magnetic  line  source.  The  slab  Is  Infinite  and  planar  with  a  thickness  2d  and  bounded 
on  both  sides  by  free  space.  The  Infinite  magnetic  field  is  applied  parallel  to  the 
slab  so  that  the  optical  axis  coincides  with  the  x  axis.  Maxwell's  equations  for  the 
situation  described  are  given  by  equations  (2-3)  and  (2-4)  with  =  0  and 
M  =  S  (x)  A  (s  -  b)  seen  that  file  direction  of  the  magnetic  line  source  Is 

chosen  perpendicular  to  the  optical  axis.  The  dielectric  properties  of  the  plasma  are 
given  by  the  relative  dielectric  tensor 

G“c  *  *  + 

p— O— O  ^0^0  —0—0 


where 


€„  =  1  -  (  ^  )  ^  (^-2) 
pa 


P  a 

! 

^MAGNETIC  LINE 

SOURCE 

Si  Mo 

b 

1 

*0 

Fig.  4el:  Geometry  of  the  source  excited  plasma  slab. 

and  a  Is  the  plasma  frequency.  The  geometry  of  the  problem  described  above  is 
illustrated  in  Fig.  4.1.  Two  cases  will  actually  be  considered^  one  where  the  source 
is  outside  the  slab  and  the  other  where  the  source  is  located  at  i  =  0  in  the  form  of 
a  slot  in  a  ground  plane. 


1.  The  Formal  Solution 

a:  The  magnetic  line  source  is  located  outside  the  slab. 

It  follows  readily  from  Maxwell's  equations  that»  for  the  configuration  in 
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Fig.  4. 1»  the  magnetic  field  will  have  only  a  y  component  and  the  electric  field  will 
have  only  x  and  z  components.  Hence,  only  F -modes  with  respect  to  the  z  direction 
will  be  of  interest.  The  symmetry  of  the  geometry  warrants  the  simultaneous 
consideration  of  the  usual  short  circuit  and  open  circuit  bisections  of  the  equivalent 
transmission  line  network.  This  is  illustrated  in  Fig.  4.  2.  The  propagation  factor 


SHORT  CIRCUIT  (sb) 
OPEN  CIRCUIT  (ob) 


Fig.  4.  2:  Equivalent  transmission  line  network  for  the  geometry  of  Fig.  4.1. 


and  characteristic  impedance  of  the  equivalent  transmission  line  in  the  region  0  <  z  <  d 
are  obtained  from  equations  (2-12c)  with  =  1  »-nd  "  6p’  Thus, 


and 


(4-3) 


(4-4) 


where  ^  is  the  wavenumber 
meters  are 


the  X  direction.  For  z>  d, 


the  corresponding  para- 

(4-5) 

(4-6) 


The  magnitude  of  the  voltage  generator  in  Fig.  4.  2  is  obtained  from  equation  (2- 22c). 
Thus,  with 


h(x,  ?  )  e  . 

vZn 

one  obtains 

V  (a,  ?  )  =  ---  6  (z  -  b). 
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The  magnetic  field  is  most  readily  expressed  in  terms  of  the  current 
function  which  satisfies  the  transmission  line  In  Fig,  4,2.  The  x  and  z  components 
of  the  electric  field  are  then  readily  derived  as  derivatives  of  These  results  can 
then  be  expressed  as  follows: 


H  f  Kz.  ?) 


.-J?* 


d?  t 


and 


E  =  -  , 


3H 


I 


aH 


E  =J-L- 
*  j  U)  e 


o 


,  0  <  ■  <  d> 


,  d  <  z  <•  , 


,  0  <  z  <  •. 


(4-7) 

(4-8a) 

(4-8b) 

(4-8c) 


The  solution  for  the  transverse  current  function  is  read!  y  obtained^  and  is  given  by 


i(*.  ?) 


.-Jkz  (b-d) 

cos  H.S 

z 

com  a.d 

■  .  (•  b). 

(4-9a) 

-jk,  (b  -  dl 

aln  a.® 

e  * 

i-k' ' ' 

B 

.  (o  b), 

(4-9b) 

•lnn^  d 


for  0  <  z  <  d,  and 


zyzTk  L  J 


for  z  >  d,  where 


(d)  = 


J - tana  d  +  — 5- 

«»«o«p  *  “«o 


H  ^ 

-J  *  cot  a  d  + 

'me  c  z  w  « 

*  orp  *0 


(•  b)  , 


(o  b) , 


(4-lla) 


(4-llb) 


and 
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T  (d)  =  << 


J—  tanned  - 

_^E__ - 

3^  H^d+k^ 


j  —  K_  cot  Kd  +  k 

K—  Z  Z  Z 


f- 


K_  cot  K_d  -  k 

**  e  Z  Z  Z 

P 
V. 

b:  The  magnetic  line  source  is  located  at  z  =  0, 


(s  b) 


(4.12a) 


(o  b) 


(4-12b) 


When  the  magnetic  line  source  is  located  at  z  =  0  in  the  form  of  a  slot  in  a 
conducting  ground  plane>  the  only  difference  in  the  equivalent  transmission  line  net¬ 
work  in  Fig.  4.2  is  that  the  voltage  source  is  located  at  z  =  0.  The  current  function 
is  then  given  by  (remembering  that  only  the  (s  b)  case  is  of  interest) 

cos  K^Z-d)  -j  k^  (z-d)l 

-E - : - j: - 


(k^  COB  K^d+  j  —  Binned) 


(4-13) 


for  0  <  z  <  d»  and 


I(z,?)  - - ^ 


-Jk^{z-d) 


K—  e  k  COB  x_d  +  Jh_  »in  x_d 
tt  P  z  z  z  z 


(4-14) 


for  d  <  z  <  »  . 

c*>  Evaliiation  of  the  free  space  field. 

In  evaluating  the  fields,  interest  will  be  centered  on  the  region  above  the  slab, 
so  that  equations  (4-10)  and  (4-14)  will  be  the  functions  in  the  integral  given  by  (4-7). 

The  treatment  of  the  fields  within  the  slab  can  be  handled  by  the  same  techniques 
presented  in  the  subsequent  outline,  although  this  will  not  be  carried  out  here.  Tamir^^ 
has  evaluated  the  fields  within  the  slab  for  the  isotropic  plasma  and  a  complete 
development  of  the  techniques  involved  can  be  found  in  his  work.  The  methods  employed 
in  the  evaluation  of  the  integral  in  (4-7)  are  well  known^*^^*  ^^and  will  be  considered 
only  briefly  here. 

To  insure  the  convergence  of  the  integral  and  the  satisfaction  of  the  radiation 

condition, it  is  required  that  Rek  >  0  and  Imk  <0.  The  integrand  has  l^ranch  points 

z  z 
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at  ?  “  ^  choice  of  branch  cuts  which  aatiaflea  the  neceaaary  conditlona  la 

c, 


Im  k^>o 
R«k^>o 


(2) 


Im  k2<  0 

R«  kj>o 


(3) 


<i) 

Rekj>o 


PATH  OF  INTEGRATION 


I"'".*"  Ml 


Fig.  4.  3:  Illustrating  the  choice  of  branch  cuts  in  the  ^  plane. 

illustrated  in  Fig.  4.  3.  The  situation  )n  Fig.  4.  3  is  depicted  for  the  case  of 
vanishingly  small  lossest  so  that  the  path  of  integration  is  clearly  defined. 

In  order  to  obtain  a  steepest  descent  approximation  for  the  fields  It  is  con¬ 
venient  to  put  the  integral  in  a  more  suitable  form.  This  is  accomplished  by  intro¬ 
ducing  a  cylindrical  coordinate  representation  and  a  change  of  variable  given  by  the 
transformation 


5  =  sin  CO  f  (4-15a) 

k  =  k  cos  (p  »  (4-15b) 

z  o 

cp  “  j  5p£  •  (4-15c) 

The  cylindrical  coordinate  system  for  the  line  source  outside  the  plasma  slab  is 
illustrated  in  Fig.  4.  4,  while  the  geometry  for  the  source  located  at  z  0  is  illustrated 
in  Fig.  4,  5.  In  the  former  case  the  change  to  cylindrical  coordinates  is  accomplished 
by  the  transformation 

X  =  r  sin  6  #  (4-16a) 

z+b-^d  -  T  cos  0  ,  (4-l6b) 

while  in  the  latter  case  the  required  transformation  is  given  by 
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z 


X  =  r  sin  0  p  (4-17n) 

z  -  d  =  r  cos  0«  (4-l7b) 

The  integrals  which  represent  the  magnetic  field  above  the  slab  can  then  be  written  as 

I  -Jk  r  cos  (cp- 0)  I 

Hy  =  -  [tt  H5f>(k^p)  -  J  e  ®  d,ipj  (4-18) 

P 

for  the  case  where  the  line  source  is  located  above  the  slab»  and 


H 


2n 


j 


-Jkr  CO*  (cD  -  6) 


Cp  CO*  k^d  CO*  «p)  +  J  *in  v) 


dcp  (4-19) 
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for  the  magnetic  line  source  located  at  z  =  0.  The  path  of  integration  in  the  cp  plane 
is  indicated  as  p«  The  transformation  given  in  equation  (4-15)  Miaps the  two  sheeted 
^  -  plane  into  a  connected  strip  in  the  op  -  plauie.  The  result  Is  indicated  in  Fig*  4.6, 
with  the  numbers  indicating  those  quadrants  In  the  C  plane  which  map  Into  the 
corresponding  strips  in  the  op  plane. 


Fig.  4.  6:  Paths  of  integration  in  the  cp  -  plane. 

The  integrals  in  equations  (4-18)  and  (4-19)  have  saddle  points  at  qp  =  0,  so 
that  a  steepest  descent  procedure  is  used  to  obtain  asymptotic  expressions  for  the 
magnetic  field.  The  steepest  descent  path  (S  DP)  is  given  by 


cos  (epy  -  0  )  cosh  ==  1. 


(4-20) 


With  the  line  source  located  outside  the  slab  the  magnetic  field  is  given  asymptotically 


'"SoT 

H  >  -  -r-^\v 

y  L 


-21) 


where 
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'  j  tan  (kd  cos  0)  - 

y  ~  » 

j  tan  (kd  cos  0)  + 


T(0)  A 


j  cot  (kd  cos  0)f  7^ 
J  cot  (kd  cos  0)  - 

’  r  H 


(s  b)  , 


(o  b). 


(4-22 a) 


(4-22b) 


In  addition  to  the  assymptotic  term  given  in  equation  (4-20)  the  field  also 
contains  residue  terms  due  to  poles  which  are  located  between  the  original  path  of 
integration  P  and  the  steepest  descent  path(S  D  F)l  A  typical  residue  term  is  given  by 


H 


yp 


(«€o 

■ifri 


- \  (<p  -  0)  ' 


(4-23) 


where  N(cp)  and  D(cp)  are  the  numerator  and  denominator  of^(cp)  respectively.  When 
equation  (4-23)  is  evaluated  the  result  is 


,  Cf  -jk  r  cos  (op  -  9) 

H  =  ^  iJL]^  o  e  °  P 

yp  ”  ^  sincp 


(4-24) 


When  the  line  source  is  located  at  z  =  0  the  asymptotic  evaluation  of  equation 
(4-19)  yields 

-j(k^r-"/4) 

H  -  uiCoe^r - - - 1  -  - - (4-25) 

'  ^Up  co8(7€^  k^d  cos  0)  +  j  sin  k^d  cos  0)  k^r, 


while  the  evaluation  for  a  particular  residue  term  yields  the  result 


H  =  - 

yp  a 


o  CSC  CD 

(JL)2 _ _£ _ 

*p  cos  (^  k^d  cos  cpp) 


-Jk^r  cos  (cpp-0) 


(4-26) 


It  is  pointed  out  that  the  radiation  field  given  by  equation  (4-25)  is  considerably  less 
complicated  than  the  corresponding  term  given  in  equation  (4-21)  for  the  case  where 
the  magnetic  line  source  is  located  outside  the  slab. 

The  expressions  obtained  above  represent  the  solutions  for  the  free  space  field. 
Before  investigating  the  properties  of  the  field  it  is  necessary  to  study  .the  transverse 
resonance  equation. 
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2«  Solution  of  the  Transverse  Resonance  Equation 

In  order  to  obtain  a  more  quantitative  picture  of  the  fields  described  in  the 
last  section  it  will  be  necessary  to  study  the  behavior  of  the  poles  of  the  various 
integrands.  These  poles  correspond  to  the  roots  of  the  transverse  resonance 
equationy  given  by 


tan  H  =  0  ,  (a  b)  , 

®P 

and 

(4-27a) 

-J  ^  cot  =  0  ,  (o  b  )  , 

aubject  to  the  auxiliary  condition 

(4- 27b) 

-1  -  •A  -  ° 

B  P  ^ 

(4-28) 

obtained  from  equations  (4-3)  and  (4-5).  For  convenience*  the  change  of  variables 

P=a,d  =  p^  +  jp^ 

q  =  .  k^d  =  +  jqj 

is  introduced.  Equations  (4-27)  and  (4-28)  then  take  the  form 

q  =  j  —  p  tan  p  ,  (a  b)  , 

*P 

(4-29a) 

q  ■  -  J  —  p  cot  p  ,  (o  b)  , 

®P 

(4- 29b) 

2  2  . 
p  -  *pq  =  0  . 

(4-30) 

By  eliminating  q  from  equations  (4-29)  and  (4-30)  one  readily  obtains  the 
equations 


1  +  —  tan^p  =  0  *  ( s  b)  * 

•p 

(4-3la) 

1  +  cot^p  =  0  *  (o  b)  • 

(4-3lb) 

It  follows  readily  from  equations  (4-31)  that  p  is  real  only  if  <0.  However*  this 
fact  will  fall  out  as  a  special  case  of  the  subsequent  development.  By  separating  the 
real  and  imaginary  parts  of  equations  (4-31)  one  readily  obtains  the  results 
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"  sln^Z  p  -  flinh^Z  p. 

-  - - - — 2  f  (  B  b)  , 

(cos  2  p^  4  coBh  2  p^) 

^  8in^2  p  -  Binh^2  p. 

- - - — 2  »  (o  b)  , 

(cob  2  p^  -  coBh  2  p^) 

and 

sin  2  p^  slnh  2  p^  =  0, 

for  both  the  (b  b)  and  (o  b)  cases.  Equation  (4-33)  permits  two  possible  casesy  namely 
either 

Pi  =  0, 

in  which  case  surface  wave  resonances  exists  or 


(4- 32a) 

(4- 32b) 

(4-33) 


P  =  -5? 

^r  2 


n  —  Of  If  2f  a  a  a 


in  which  case  leaky  wave  resonances  exist. 

In  the  surface  wave  case*  substitution  of  p^  =  0  in  equations  (4-32)  yields 


e 


P 


*  2 

-  tan 

-  cot  Py 


»  (s  b), 
•  (o  b). 


(4- 34a) 
(4- 34b) 


It  is  seen  that  the  existence  of  surface  waves  is  possible  only  when  <0,  that  iSf 
when  the  plasma  is  opaque.  For  the  case  where  p^  ^^^12,  equations  (4-32)  yield 
the  result 


e 


P 


-  tanh  p 


i 


(4-35) 


for  both  bisectionSf  with  the  stipulation  that  n  is  even  for  the  (s  b)  case  and  odd  for 
the  (o  b)  case.  The  above  results  can  be  summarized  by  sketching  tlie  loci  of  the  poles 
in  the  p-plane  with  frequency  as  a  parameter.  The  sketch  is  shown  in  Fig.  4.  7.  It 
is  seen  that  the  surface  wave  pole  loci  lie  along  the  p^  axis  and  the  surface  wave 
poles  become  complex  at  the  plasma  frequency.  It  is  also  noted  that  the  surface  wave 
poles  are  infinite  in  number  for  any  frequency  below  the  plasma  frequency. 
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Fig*  4«  7:  Pole  loci  in  the  p-plane«  Arrows  indicate  increasing  frequency. 

The  roots  in  the  q-plane  can  be  obtained  from  equations  (4*29).  ThuSf  in  the 
case  of  the  surface  waves  with  p  =  p^  and  q  =  jq^  one  obtains 


-p,  cot  , 
But  from  equation  (4-30)  one  can  write 


P, tan  p^ 


*p**l 


(a  b)  , 
(o  b)  . 


(4-36a) 
(4- 36b) 


(4-37) 


and  hence  the  result 
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'  -Py  cot  P,  .  (■  b)  ,  (4- 38a) 

qi  =  ■ 

[  Pytanp^  ,  (ob)  ,  (4- 38b) 

is  obtaineH,  Equations  (4-38)  can  be  used  to  obtain  the  surface  wave  loci  shown  in 
Fig.  4.8. 


q=Jq| 


n  =  0 

tb 


PROPER 

RESONANCES 


IMPROPER 

RESONANCES 


Fig.  4.  8:  Surface  wave  pole  loci  In  the  q- plane. 

In  order  to  obtain  the  loci  of  the  leaky  wave  poles  in  the  q-plane^  equation 
(4-35)  and  the  fact  that  p^  ==  ^/2  are  substituted  into  equations  (4-29).  This  last 
operation  yields  q  as  a  complex  function  of  frequency  (or  f^)g  and  after  separating 
the  real  and  imaginary  parts  and  eliminating  the  frequency  parameter^  one  obtains 
for  the  equation  of  the  leaky  wa\re  loci  in  the  q-plane. 

2qr  .1  2q 

q  =  -  _  coth  —  ,  n  =  1|  Z,  3p . 

1  nn  nn 

The  range  of  values  on  q^  is  given  by 


(4-39) 
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q  =  dfc  •  #  ifor  uo  -  w  or  f  s  0, 
r  P  P 

and 

i  ^  #  for  a  =  •  or  Cp  =  L 

For  the  special  case  where  n  ^0,  equations  (4-29a)  and  (4«-36)  yield 

qi  =  -  Pi  coth  p^,  (4-40) 

in  which  case  takes  on  values  from  -1  to  -«  with  increasing  frequency.  It  is  pointed 
out  that  a  pole  given  by  equation  (4-40)  will  not  contribute  to  the  field.  This  fact  will 
become  apparent  when 


Fig.  4.  9:  Lieaky  wave  loci  in  the  q-plane. 

the  pole  loci  are  mapped  into  the  op  plane.  The  loci  of  the  leaky  wave  poles  are 
sketched  in  Fig.  4.  9  with  the  arrows  indicating  increasing  frequency. 

The  mapping  of  die  pole  loci  into  the  ^  -  plane  is  accomplished  by  means  of  the 
transformation  given  by  equation  (4-15b).  Thus^  one  obtains 

Q  *-k  d  =  -  k  d  cos  op  I  (4-41) 

^  a  o 


with 


«p  =  ♦  J«pi  • 
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E^uatloD  (4-41)  li  readily  ■•pirated  into  the  real  and  Imaginary  parts  to  yield 


qr  =  -  k^d  coe  cosh  tp^  ,  (4-42a) 

=  k^d  sin  tp^  sinh  cpj  .  (4-42b) 


The  eurface  arave  loci  In  the  q-plane  were  all  located  on  the  Imaginary  axle 
and  hence  It  ie  recognised  from  equation  (4-42a)  that  these  loci  will  lie  on  the  line 
(py  -  ^ IZ  In  the  (p  -  plme  .  Equation  (4-42b)  can  then  be  written  as 


■inh 


«  -  *** 
«Pl  -T^ 


(4-43) 


Equation  (4-43)  can  bo  wrltteng  with  tho  aid  of  equation  (4-30)g  ao 

Z 

I 

•inh  q 


(4-44) 


where  c  ie  the  apeed  of  light.  *rhen,  considering  ttie  n  =  0  resonance  it  is  seen  that  as 
iS  approaches  serog  p^  approaches  ^/Z  and  equation  (4-44)  yields 


.  ,  ^  _  TTC 

"Pi  ~  •  liTT  • 

P 

Alsog  as  u)  approaches  p^  approaches  zero  and  approaches  -Ig  so  that  equation 
(4-43)  yields 


•inh 


"  ■  u»pT  • 


It  is  noted  that  a  surface  wave  pole  for  the  n  -  0  case  cannot  contribute  a  residue  term 
to  the  fields.  By  reasoning  similar  to  the  above  one  can  obtain  the  behavior  of  the 
higher  surface  wave  loci  in  the  op  -  plane.  Thusg  for  the  next  three  resonancesg  the 
following  results  are  obtained: 


n  =  Ig 
at  (1)  -  0 


at  ttj  s  dip 


•inh  cp^  = 


0 


TTC 


(proper  resonance) 
(improper  resonance)* 


•inh  cp^  B  « 
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n  -  2| 
at  *  0 


•inh 


«Pi 


at  uu  -  (JD 

n  =  3 
at  oy  =  0 


ainh  i 


TT  C 

2  (ju  d 

P 

3ttc 

ZF^ 

ainh  ±  • 

TT  C 

(a  d 
P 

2lic 


(proper  reaonance) 
(improper  reaonance)^ 


(proper  reaonance) 
(improper  reaonance)t 


at  tt)  =  (a^ 


ainh  op^  =  ±  «  • 

The  above  information  for  aurface  wave  polea  ia  aummarised  in  Fig.  4.10. 
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Fig.  4.  10:  Surface  wave  pole  loci  in  the  op  -  plane 
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The  location  of  the  leaky  wave  loci  in  the  cp  *  plane  is  a  more  complicated 
matter.  In  this  case 


p  =  i  i  j  tanh  ^  , 


and 


=  ;  -£IL_.  j_L 


Equations  (4*42)  can  then  be  written 


(4*45) 


(4-46) 


k^d  cos  ep  ^  cosh  CP  ^  =  i 


nn 


k^d  sin  cpy  slnh  cp^  =  -  -  ^ —  tanh*^  » 


(4-47a) 


(4*47b) 


1  -1 _ 

From  equation  (4- 47b)  It  can  be  concluded  that»  since  *—  tanh  >  0#  leaky  wave 

jr,  ’’ 

poles  will  map  into  the  strips  0  <cp^  <  rr  *  cp^  <  0  and  *n<  cp^  <  0^  cp^  >  0.  Hence, 
contributions  to  the  fields  will  be  due  only  to  leaky  wave  poles  located  in  the  region 
0  <  cp^  <  ^/2,  cp^  <0.  Also,  from  equation  (4-47b),  it  is  seen  that  the  poles  are 
located  symmetrically  about  the  origin  of  the  cp  -  plane*  In  subsequent  calculations 
attention  will  be  limited  to  the  strip  0  <cp^  <  ^/2,  cp|  <  0. 

Next,  limiting  values  will  be  calculated*  Equations  (4*47)  can  be  written  as 


cos  cp^  cosh  cp^  - 


nnc 


(4* 48a) 


sin  cp  ^  sinh  cp  ^ 


tanh"^  Jl  -  (*p/ti))^ 

'7 


(4*46b) 


<1 J 


Also  from  equations  (4*47)  one  obtains 


,  UTT 


\L"r 


coth  -^tan  cp^  tanh 


'Pi]' 


(4-49) 


Now,  when  q^  approaches  *^^/2  it  follows  from  (4*49)  that 
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tan  tanh  Qp^  s  «  , 


and  hence  cp  -  2*  corresponding  value  of  cp^i  obtained  from  equation  (4-48b)  is 

=  0.  These  limiting  values  correspond  to  m  -  When  u)  =  u)p  the  limiting  values, 
obtained  from  equations  {4-48b)  and  (4-49)  arc  given  by  =0,  . 

Further  properties  of  the  roots  can  be  obtained  from  equation  (4- 48a).  Note 
that  at  a  given  finite  frequency,  as  n  is  made  to  increase,  cosh  must  increase  (since 
0  <  cos  (p  <  1  in  the  region  of  interest).  This  means  that  the  poles  corresponding  to 
large  n  at  this  fixed  frequency  must  be  in  the  region,  cp^  small  and  cp.  very  negative. 
Consequently,  it  follows  (from  the  shape  of  the  steepest  descent  curve  in  Fig.  4.  6) 
that  a  finite  number  of  leaky  wave  poles  will  contribute  to  the  field  for  a  given  angle 
of  observation  0  <  0  <  tt/2. 

The  above  information  yields  a  qualitative  picture  of  the  leaky  wave  pole 
loci  in  the  qj  -  plane.  It  is  possible  to  obtain  a  more  exact  picture  by  graphical  calcula¬ 
tion.  Thus  equations  (4-47)  and  (4-49)  can  be  used  to  obtain  an  explicit  equation  for 
the  pole  loci,  namely 


(4-50) 


where  \  is  the  plasma  wavelength.  Equation  (4-50)  can  be  solved  graphically  to  yield 
p 

quantitative  reaults  for  the  poles  In  the  tp  -  plane.  The  qualitative  Information  obtained 
above,  together  with  a  knowledge  of  the  quantitative  behavior  of  the  pole  loci  obtained 
from  equation  (4-50),  shows  that  the  leaky  wave  loci  In  the  cp  -  plane  are  as  Illustrated 
In  Fig.  4-11  a  and  b.  The  steepest  descent  path  for  6  ="12  is  shown  dotted. 


Fig.  4.Ua: 
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3.  Some  Properties  of  the  Surface  Waves 

One  of  the  unexpected  features  of  the  surface  wave  resonances  is  the  fact 

that  an  infinite  number  of  these  exist  and  give  rise  to  residue  terms  at  any  frequency 

below  the  plasma  frequency*  Actually  this  result  is  not  so  strange  when  one  considers 

the  dispersion  relations  given  by  equations  (4-3)  and  (4-5).  These  are  plotted  in 

Fig*  4*  12»  on  the  same  set  of  axes^  for  u)  <  U)  •  It  is  seen  from  an  examination  of  the 

curves  in  Fig*.  4*  12  that  an  infinite  number  o/’imaglnary  values  of  k  will  satisfy  both 

z 

curves  for  real  values  of  ^  and  Another  point  of  interest  is  the  fact  that  the 
vector  drawn  from  the  origin  to  a  point  on  the  plasma  dispersion  curve  (i*  e*  the 
wavenumber)  tends  to  become  collinear  with  the  line  k  =-—1—  5  for  higher  surface 

*  l‘p| 

wave  resonances  (higher  values  of  ^)*  It  Is  also  noted  that  for  the  higher  resonances 
the  values  of  g  and  larger  and  hence  the  electrical  length  of  the  slab  increases 

with  n« 

The  dispersion  curves  for  the  individual  surface  waves  can  be  easily  obtained 
from  the  results  which  have  already  been  derived*  Hence  *  using  the  fact  that 
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Fig.  4.  12:  Plot  of  dispersion  relations  for  free  space  and  for  the  plasma  when  m< 
and  introducing  the  results  given  in  equations  (4-34)  one  obtains 


-  sec  p^  »  (s  b)  g 

CSC  Py  »  (o  b) 


Nowg  introducing  equations  (4-38)  into  equation  (4-5)  one  obtains 


(4-51a) 
(4- 51b) 


Equations  (4-51)  and  (4-52)  are  the  parametric  equations 


,  (s  b)  ,  (4-52a) 

g  (o  b).  (4-52b) 

of  the  surface  wave  dispersion 


curves*  The  range  of  values  on  the  parameter  p^  can  be  obtained  from  Fig,  4.  7,  with 


the  result  that  the  dispersion  curves  can  be  sketched  as  shown  in  Fig*  4.13. 


35 


The  elopes  of  the  curv^es  in  Fig.  4*13  are  all  negative.  Since  the  group  velocity  is 
related  to  the  negative  of  the  elopes  of  these  curve s^  the  phase  and  group  velocities 
are  both  positive  and  the  waves  are  all  forward  in  the  x  direction. 

The  spatial  variation  of  the  surface  waves  is  sinusoidal  within  the  slab  and 
exponentially  decaying  outside  the  slab.  Ift  in  the  source  free  case#  is  the  amplitude 
of  the  current  function  at  z  =  d  the  field  distribution  can  be  written  as  follows: 


s  b  case 
( a)  0  <  z  <  d 
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(4-53a) 
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(b)  z  >  d 


vth;  ^“^o'p  ‘=°'  »‘  z‘* 


^/TTT 


(4-53b) 
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E  =.A__L. 


^2n 


U)  e. 
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•InR  d 

X 


(b)  z>  d 


E  = 
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cos  K 


jn^  ^®*o'p 


d  -Ik  |(z  -  d) 

«  .-j?*,  I  »r 


I.  -  k  (z-d) 

H  =— ^  •  '  *1  ,  (4-54b) 

^  yziT 

E  =-!s_  -5.  .-)?». -N 

‘  /IT  “‘o 

One  of  the  interesting  consequences  of  the  anisotropy  in  the  present  problem 
is  that  the  s  component  of  the  electric  field  is  continuous  over  the  plasma- free  space 
interface.  As  a  result  power  travels  in  the  same  direction  both  inside  and  outside  the 
plasma.  This  result  follows  from  the  fact  that  the  only  real  component  of  the  Poynting 
vector  is  s  -  ^x^y  *  interesting  to  point  out  that  in  the  case  of  the  isotropic 
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plasma  slab#  which  has  been  thoroughly  studied  by  Tamir^\  power  travels  in  opposite 
directions  in  the  plasma  and  free  space  regions.  Againi  this  result  occurs  because 
of  the  discontinuity  in  over  the  interface. 

One  more  point  of  interest  should  be  mentioned  in  connection  with  power.  If 
the  fields  within  the  slab  are  interpreted  as  being  composed  of  plane  waves  which  are 
reflected  back  and  forth  within  the  slab,  then  the  wavenumber  for  a  particular  plane 
wave  is  given  by  the  vector  which  is  constructed  in  Fig,  4.12.  The  average  power 
flow  vector  P  =  (5.  *  U  ^  directed  normal  to  the  dispersion  curve®*  It  is 

seen  from  Fig.  4.12  that  the  angle  between  the  wavenumber  and  P  is  less  than  90”.  It 
is  pointed  out  that  the  net  power  transfer  is  in  the  x  direction  even  though  Phas  a 
component  in  the  negative  z  direction.  This  follows  from  the  fact  that  when  the  plane 
wave  is  reflected  at  the  plasma  interface  the  resulting  Poynting  vector  will  have  a 
component  in  the  positive  z  direction.  Hence,  the  Poynting  vectors  due  to  the  two 
plane  waves  combine  to  give  a  net  component  in  the  x*direction. 


4.  Some  Properties  of  the  Fields  Above  the  Slab 

Some  of  the  properties  of  the  fields  due  to  the  residue  terms,  and  their 
relationship  to  the  geometry  of  the  structure  will  now  be  considered.  Since  the  general 
properties  of  the  surface  and  leaky  wave  fields  are  essentially  independent  of  the 
source  location  the  discussion  will  be  limited  to  the  case  where  the  line  source  is 
located  at  z  =  0.  Also,  since  the  characteristics  of  such  waves  have  been  discussed 
in  great  detail  elsewhere,  ^Hhe  presentation  here  will  be  brief. 

From  equation  (4-26)  it  is  seen  that  the  residue  fields  vary  as 


where 


-Jk  r  cos  {»  -0) 
f  (r.  e)  =  e  ^ 


Vp  -  -T  +  J  % 


(4-55) 


for  the  surface  waves,  and 


'Pp  =  'Prp  i'Plp 

for  the  leaky  waves.  Now,  considering  the  surface  waves  separately,  it  is  seen  that 
the  constant  phase  surfaces  are  perpendicular  to  the  x-axls  and  given  by 


r  cos  -  0)  =  constant, 

while  the  constant  amplitude  surfaces  are  parallel  to  the  x-axis  and  given  by 


r  sin  (j  -  0)  =  constant. 
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It  follows  that  if  a  surface  wave  is  observed  along  a  plane  of  constant  phase  the 
amplitude  decreases  exponentially  with  the  distance  from  the  interface.  These  remark?* 
are  illustrated  in  Fig.  4. 14.  The  angle  0^  is  the  angle  of  observation  at  which  the 


Fig.  4.14:  Constant  phase  and  amplitude  surfaces  for  the  surface  wave  field  above 
the  slab. 


residue  field  for  a  particular  surface  wave  pole  is  first  observedf  i,  e.  t  the  angle 
for  which  the  steepest  descent  ciurve  intercepts  the  pole. 

Remarks  which  are  similar  to  the  above  can  also  be  made  for  the  leaky  waves. 
In  this  case  ep^  <  6  and  the  constant  phase  surfaces  are  given  by 


r  cos  (0  -  ffi  )  =  constantt 
Trp 

while  the  constant  amplitude  surfaces  are  given  by 

r  sin  (0  -  n  )  =  constant. 

^rp 

It  follows  from  these  results  and  equation  (4*55)  that  the  leaky  wave  decreases 
exponentially  with  distance  from  the  wave  vector  k^  which  is  indicated  in  Fig.  4. 15. 
This  last  result  implies  growth  in  the  direction  perpendicular  to  the  slab  and  decay 
in  the  direction  parallel  to  the  slab.  These  remarks  are  illustrated  in  Fig.  4.15. 

One  more  fact  is  of  interest  and  this  follows  from  the  expression  for  the  phase  velocity 
for  a  surface  or  leaky  wave 


-  c 

'"p  ^  cot  tS-epyp)  coilifjp* 
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Fig.  4. 15:  Constant  phase  and  amplitude  surfaces  for  the  leaky  wave  field  above  the 
slab. 

It  is  seen  that  the  surface  waves  are  slow  waves,  while  the  leaky  waves  are  fast  waves. 


« 


40 


REFERENCES 


L  R.  E.  Collin,  Field  Theory  of  Guided  Waves,  McGraw  Hill  Book  Company,  Inc. , 
pp.  97-116,  New  York,  1960. 

2.  R.  E,  Collin,  A  Simple  Artificial  Anieotropic  Dielectric  Medium”,  I.  R.  E.  Trans., 

vol.  MTT-6,  pp.  206-209;  April,  1958.  - 

3.  G.  Joo8,  Theoretical  Physics,  2d  ed. ,  Blackie  and  Sons,  Ltd.,  Glasgow,  1951. 

4.  G.  Estrin,  "The  Effect  of  Anisotropy  in  a  Three-Dimensional  Array  of  Conducting 

Disks”,  Proc.  I.R.E.,  vol.  39,  pp.  821-826;  July,  1951. 

5.  G.  Estrin,  ”The  Effective  Permeability  of  an  Array  of  Thin  Conducting  Disks”, 
Journal  of  Appl.  Phys. ,  vol.  21,  pp  667-670;  July,  1950. 

6.  E.  Arbel  and  L.  B.  Felson,  ”On  Electromagnetic  Green’s  Functions  for  Uniaxially 
Anisotropic  Regions”,  Electrophysics  Dept. ,  Polytechnic  Institute  of  Brooklyn, 
Report  PIBMRl-985-61;  Feb.  ,  1962. 

7.  E.  Arbel,  "Radiation  from  a  Point  Source  in  an  Anisotropic  Medium”,  Microwave 

Res.  Inst.,  Polytechnic  Institute  of  Brooklyn,  Report  PIBMRI-861-60;  Nov.,  I960. 

8.  L.  B.  Felson,  "Radiation  from  a  Uniaxially  Anisotropic  Plasma  Half  Space”, 
Electrophysics  Dept. ,  Polytechnic  Institute  of  Brooklyn,  Report  PIBMRI-1058-62; 
August,  1962. 

9.  L.  B.  Felson  and  N.  Marcuvitz,  "Modal  Analysis  and  Synthesis  of  Electromagnetic 

Fields”,  Polytechnic  Institute  of  Brooklyn,  Microwave  Res.  Inst.  ,  Reports 
R-446-55  (a,b),  R-726-59,  R-776-59. 

10.  S,  Barone,  "Leaky  Wave  Contributions  to  the  Field  of  a  Line  Source  Above  a 
Dielectric  Slab”,  Polytechnic  Institute  of  Brooklyn,  Microwave  Res.  Inst.  , 

Report  PIBMRI  532-56;  Nov.  1956. 

11.  T,  Tamir,  "The  Electromagnetic  Field  of  a  Source  Excited  Plasma  Slab”, 
Polytechnic  Institute  of  Brooklyn,  Microwave  Res.  Inst.,  Report  PIBMRI-932-61; 
Sept.  1961,* 

12.  T,  Tamir,  "Leaky  Wave  Contributions  to  the  Field  of  an  Electric  Line  Source 
Above  a  Plasma  Slab”,  Polytechnic  Institute  of  Brooklyn,  Microwave  Res.  Inst,  , 
Report  PIBMRI-845-60;  Nov.  1960. 

13.  T.  Tamir,  "The  Influence  of  Complex  Waves  on  the  Radiation  Field  of  a  Slot- 
Excited  Plasma  Layer",  Polytechnic  Institute  of  Brooklyn,  Microwave  Res.  Inst. , 
Report  PIBMRI-959-61;  Nov.  1961. 

14.  B,  Friedman,  Principles  and  Techniques  of  Applied  Mathematics,  J.  Wiley  and 
Sons,  Inc,,  New  York,  1957. 

15.  L.  D,  Landau  and  E.  M.  Lifschltz,  Electrodynamics  of  Continuous  Media  , 
Pergamon  Press,  New  York,  I960. 

16.  J,  G.  Linhart,  Plasma  Physics,  Intersclence  Publishers,  Inc,,  New  York,  I960. 


t 


41 

17.  K.’G.  Buddent  The  Wave-Guide  Mode  Theory  of  Wave  Propagation#  Prentice - 
Halit  Inc. ,  New  Jersey,  1961. 

18.  A.  Erdelyi,  Asymptotic  Expansions,  Dover  Publications,  Inc. ,  New  York;  1956. 

19.  N.  Marcuvitz  and  J.  Schinger,  *'On  the  Representation  of  the  Electric  and 
Magnetic  Fields  Produced  by  Currents  and  Discontinuities  in  Waveguides'*, 

J.  Appl.  Phy.,  22,  pp  806-819,  195L 

20.  K.  G.  Budden,  Radio  Waves  in  the  Ionosphere,  Cambridge  University  Press,  1961. 


